Dynamical Crystallization in the Dipole Blockade of Ultracold Atoms 
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We describe a method for controlling many-body states in extended ensembles of Rydberg atoms, 
forming crystalline structures during laser excitation of a frozen atomic gas. Specifically, we predict 
the existence of an excitation number staircase in laser excitation of atomic ensembles into Rydberg 
states. Each step corresponds to a crystalline state with a well-defined of regularly spaced Rydberg 
atoms. We show that such states can be selectively excited by chirped laser pulses. Finally, we 
demonstarte that, sing quantum state transfer from atoms to light, such crystals can be used to 
create crystalline photonic states and can be probed via photon correlation measurements. 

PACS numbers: 32.80.Rm, 03.67.-a, 42.50.Gy 



When cold atoms are excited into high-lying Rydberg 
states, the resulting interactions give rise to energies 
that exceed the translational energy by many orders of 
magnitude and strongly modify the excitation dynamics. 
In very small ensembles confined to a few micrometers, 
a single Rydberg atom can entirely block any further 
excitation. This, recently observed, "dipole blockade" 
[l], 0, IE Ej], enables the production of highly entangled 
collective states with potential applications for fast quan- 
tum information processing [E H^Sj], and as single- atom 
as well as single-photon sources [2] . 

A number of theoretical and experimental studies have 
addressed the evolution of Rydberg states in l arg e atomic 
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strong interactions are critical for many of the predicted 
and observed phenomena, they also cause a very com- 
plex excitation dynamics, having, thus far, prevented 
a controlled quantum state preparation in larg e ensem- 
bles. For example, experiments [7L lil @, [ToL 111! Il2j have 
probed interesting relaxational excitation dynamics to- 
wards equilibrium states whose releva nt p hysics is well 
captured by mean- field descriptions @, EE EH . 

In this work we describe a method that overcomes this 
obstacle and permits coherent manipulation of strongly 
correlated, many-body states of large ensembles. The ap- 
proach is based on a dynamical crystallization of local- 
ized, regularly-spaced collective excitations. The crys- 
talline states can be created with chirped laser pulses, 
and have a widely tunable lattice spacing. By extend- 
ing the dipole blockade mechanism [2] to the regime of 
multiple Rydberg atoms, such excitation crystals provide 
a suitable initial state to realize quantum random walks 
EE 0] • We show that the crystalline correlations can be 
mapped onto a propagating light beam and be detected 
via photon correlation measurements as a pulse train of 
localized single photons. 

The key idea is illustrated in fig. [TJ For repulsive van 
der Waals interactions and for a given Rydberg atom 
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FIG. 1: Adiabatic crystal state preparation, (a) A cold 
atomic gas is illuminated by a chirped laser pulse, whose fre- 
quency changes from the blue (A < 0) to the red (A > 0) side 
of the driven Rydberg transition, (b) Schematics of the many- 
body energy spectrum as a function of the laser detuning A 
and for a Rabi frequency Q = 0. Shown are the lowest-energy 
states \n) for a given Rydberg atom number n. These states 
correspond to Rydberg atom (green spheres) lattices sketched 
in (c). Their energy crossings give rise to an excitation num- 
ber staircase shown in (d) for a chain of N — 31 atoms (blue 
curve) and a cigar-shaped condensate (red curve) with the 
same length s\\ and 500 atoms. The insets show the Rydberg 
atom densities resulting in the latter at the detunings marked 
by the arrows. Panel (e) shows the ground state energy gap 
for the chain (blue curve) and the corresponding gap distri- 
bution of the condensate (red shading) at O = 10 3 • C/sn. 



number n, the states with the lowest interaction energy 
maximize the spacing between Rydberg atoms, which, 
hence, arrange on regular crystals (fig. []})). If initially all 
atoms are in their ground states, for large negative detun- 
ings A of the laser field, the initial state coincides with 
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the many-body ground state in the rotating frame of ref- 
erence. If A is adiabatically increased to positive values, 
it effectively lowers the energy levels of excited many- 
body states (fig. []J). They ultimately cross at critical 
detunings A™ , l5 marking discrete jumps in the Rydberg 
atom number from nton+1 (fig.UJi). As a result, larger 
and larger Rydberg atom crystals are consecutively pop- 
ulated during evolution of the chirped pulse. In essence, 
the detuning A acts as a control parameter that decreases 
the energy difference between adjacent number states |n) 
and \n + 1), allowing further Rydberg atoms to "enter" 
the excitation volume. This interaction blockade is sim- 
ilar to a recently observed effect of interacting atoms in 
an optical lattice [21]. Moreover, the resulting "dipole- 
blockade staircase" bears resemblance to the Coulomb 
blockade staircase observed in nano-scale solid-state de- 
vices [22[, where the number of electrons increases as 
a function of external lead- volt age. In contrast to such 
solid state devices, the present system permits dynamical 
manipulation of strongly interacting particles in a highly 
coherent manner. 

The Rydberg gas is described as an ensemble of two- 
level systems at fixed positions r$, each possessing a 
ground state \gi) and an excited state |e$) [Hf, laser- 
coupled with a Rabi frequency Q(t). Within the frozen 
gas limit, the excitation dynamics is described by the 
Hamiltonian 
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where <r^ 

describes the repulsive (C > 0) van der Waals interaction 
between Rydberg atom pairs at a distance . 

We first consider the simplest case of atoms in a one- 
dimensional optical lattice [23|, [24| , modeled as a chain of 
length si I , that is composed of N atoms equally separated 
by a = S||/(iV — 1). At zero laser intensity, the general 
structure of the resulting many-body energy spectrum 
can be easily understood (see fig. [Jj). The 2^ many- 
body states are composed of states |n, fc), where 

k labels the spatial configuration of n excited Rydberg 
atoms. Their partially degenerate interaction energies 

(k) 

€n are shifted by — nA. Within each n-manifold, the 
state |n,0) = \n) with the lowest energy E^ = E n = 
— nA maximizes the Rydberg atom spacing at Zij = 
s\\/(n — 1) = a n , forming a Rydberg atom crystal with 
lattice spacing a n . For small excitation fractions (a <C 
a n ), the corresponding energies 
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Tuning A from negative to increasingly positive values, 
thus, induces stepwise changes in the many-body ground- 
state character from the initial state |0) to regular chains 
of n excited Rydberg atoms (|n)). Equation (pQ) yields a 
well-defined large- N limit, where the detuning 
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determines the fraction v = n/N of excited atoms and 
the Rydberg atom lattice spacing a n = s\\/(n— 1). 

Transitions between the crystal states are induced by 
the excitation laser, which couples adjacent n-manifolds, 
such that the degeneracies at A™ +1 are lifted and re- 
placed by avoided crossings of separation . Adiabatic 
preparation of a given crystal state is possible provided 
that the pulse-evolution is slow compared to l/<$™ +1 . 

At the first three crossings adjacent lattice states are 
directly coupled by the laser field (figEH-c), such that the 
energy gap S™ +1 ~ ^. For n > 3, however, there is no 
direct laser-coupling between adjacent levels (fig.[2ji). In- 
stead, the gap arises from (2n — 3)-photon transitions via 
off-resonant intermediate states detuned by ~ C/a^ +1 . 

If at the time t of a given n — > (n + 1) transition, 
uo n = VL(t) / (C a~+i) <C 1, the intermediate states, that 
couple \n) and \n + 1) are far off-resonant and act as vir- 
tual levels for resonant multi-photon transitions. For a 
given Rydberg atom number n these transitions involve 
2n — 2 intermediate states such that the gap scales as 

K+i ~ ^n n " 3 ( n /( N ~ i)) 6 c / a6 This exponential drop 
with increasing n, prevents the preparation of large crys- 
tals. However, at larger Rabi frequencies, uj n > 1, power 
broadening exceeds the intermediate state detunings. In 
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FIG. 2: (a)-(c) Schematics of single-photon coupling between 
adjacent Rydberg lattice states \n) and \n+ 1) for n < 3. (d) 
Schematics of the three-photon coupling for n — 3. The inter- 
mediate states are off-resonant by ~ C/a Q n+1 . (e) Resulting 
energy gap for different N and n as a function of Q. The fitted 
power- laws demonstrate the transition from multi-photon to 
to sequential one-photon coupling, discussed in the text. 



yield successive level crossings between regular n- and 
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FIG. 3: (a) Pulse envelope and chirp of the applied laser 
field, that excites a cigar-shaped Rb condensate to 65s states. 
(b)-(e) Resulting time dependence of the probabilities P n to 
excite n atoms in a gas of 300 atoms and with a length of 
s\\ — 15/xm (b), 22.5/xm (c), 35/mi (d), 45/i (e). The insets 
show the resulting final Rydberg atom densities. 



this regime, consecutive crystal states are coupled by se- 
quential one-photon transitions, which results in a linear 
^-scaling of 5™ +1 . 

This is confirmed by our numerical calculations, 
demonstrating the transition from multi- to single- 
photon coupling at Q ~ C/a^ (see fig. [2^). Beyond this 
critical Rabi frequency the gap scales as ~ z/fi, en- 
suring a finite gap in the large- n limit. Since the lattice 
spacing a n ~ \j(yp) is also independent of sn, arbitrarily 
large Rydberg crystals should indeed be producible with 
finite peak Rabi frequencies ft > C/a^. 

We now focus on the excitation dynamics in 
frozen gases, of randomly spaced atoms. To be 
specific, we consider a Bose-Einstein condensate, 
whose initial state Yli 0( r i) ' ®i \9i) ls described by 
an axially symmetric Thomas-Fermi profile 0(r) oc 

Observables are now de- 
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termined by quantum average over the TV-body spatial 
wave function and evaluated through Monte- Carlo inte- 
gration. The fully correlated iV-body wave function is 
obtained via a Hilbert space truncation approach that 
exploits the excitation blockade of closely spaced atoms 

a 

The simple chain model, developed in the previous 
sections, provides a qualitative picture for the excita- 
tion dynamics in disordered systems (see fig. [TJl-e). An 
important difference is that a disordered gas accommo- 



dates Rydberg atom configurations with continuously 
distributed interaction energies, leading to unavoidable 
non-adiabatic transitions within a given n-manifold at 
the end of the laser pulse. The resulting crystals are 
thus composed of localized collective excitations. 

To demonstrate the experimental feasibility of the pro- 
posed scheme, we have simulated the dynamics of a cigar- 
shaped Rubidium condensate excited to 65s Rydberg 
states [lH with realistic laser pulses (il. fl2j| (fig. |3^). For 
s\\ = 15/im, two excitations, localized at the edges of the 
cloud, are excited. Increasing su gives smaller critical 
detunings A™ +1 [see equation (pQ)], such that consecu- 
tively larger lattices are produced (fig. 03>|3k). Further 
increasing the cloud length we find a linear increase 
of the correlation length up to our largest systems size of 
70/im. The depicted time evolution of the Rydberg atom 
number populations P n demonstrates that — despite the 
disordered atom positions — the gas is transfered to a 
state with precisely n excitations, arranged on perfectly 
filled Rydberg crystals. 

The described mechanism equivalently applies to 
higher dimensions and different geometries. Exemplary 
we show in fig. [H the excitation-number staircase of a 
pancake-shaped condensate (s± ^> sn) along with some 
Rydberg atom densities of the many-body ground state. 
The depicted densities are obtained via a Monte Carlo av- 
erage over randomly sampled atom positions, from which 
we obtain the lowest energy excited state configuration 
for a given n. To visualize the angular correlations, each 
individual configuration is rotated such that one exci- 
tation is centered at the bottom of the image and the 
horizontal component of the Rydberg atoms' center of 
mass vanishes. Analogously to the one-dimensional case, 
Rydberg atoms are excited on regular lattice structures, 
which undergo a structural transition to a single-atom 
centered hexagon, with increasing A or s±. Increasing 
the system size s± : thus, yields radial growth of large 
triangular crystals, with the lattice constant controllable 
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FIG. 4: Dipole blockade staircase in a pancake-shaped con- 
densate. The insets show Rydberg atom densities of the 
many-body groundstate at the indicated plateau values of the 
laser detuning A. 
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FIG. 5: Light mapping of a three- Rydberg atom crystal 
with S|| = 30/xm and s± = 4/rni, embedded in a conden- 
sate with a longitudinal and transverse size of 50/xm and 
10/im, respectively, (a) Pulse propagation dynamics, as de- 
scribed by the expectation value of the photon density opera- 
tor p p h oc E^E. The outer- left blue curve shows the time evo- 
lution of the control Rabi frequency Q f , while the right green 
area gives the resulting photon signal, detectable outside the 
atom cloud, (b) Corresponding two-photon correlation func- 
tion — (pph(£i)p p h(£2)). Multiplication by the squared 
speed of light, c, yields a temporal density and ensures the 
correct normalization c 2 f g^dt\dt2 = n(n — 1). 



by the detuning A. 

Single-atom decay processes, such as black body radia- 
tion induced transitions and spontaneous decay typically 
occur on very long timescales of several 100/i [28[. In 
the limit of small final decay probabilities, we can esti- 
mate its importance by following the coherent dynamics 
and determining the decay probability from the time- 
dependent number of excited atoms. As shown in Figj3] 
the non-fidelity due to such processes is as small ~ 0.02, 
even for the largest considered Rydberg chain. As this 
values increases linearly with n, the preparation of rather 
large chains of up to 25 atoms has a decay-limited fi- 
delity of T « 0.9. Cryogenic environments [29[ reduce 
the value of 1 — T by a factor of ~ 3 for the present 
parameters [28], thus, permit high fidelity preparation of 
considerably larger crystals. A less controllable decoher- 
ence source arises from atomic motion induced by the 
Rydberg-Rydberg atom interaction itself. We estimate 
its effect by the classical displacement Ar ~ CV -7 M -1 r 2 
of an atom pair, initially separated by r, where M is 
the mass of the atoms. For the parameters of fig. [3] 
(r > 10/xm) and r = 4.0/is this gives a change in atom 
distances by less than a fraction of 10 -3 , leaving the ex- 
citation dynamics practically unaffected. 

A number of techniques can be used to experimentally 
probe the Rydberg crystals. Absorption imaging of the 
groundstate atoms in optical lattices [26[ could allow to 
detect excited atoms as defect spots at the corresponding 
lattice sites. Electric field ionization of Rydberg atoms 
combined with position sensitive MCP detectors and ion 
optics magnification [27[ could allow to directly image 
the predicted /im-scale structures. Already measuring 
the final number distribution of Rydberg atoms [9[ would 
indirectly verify the production of the discussed crystals 



states, since they stand out due to a single number n of 
excited atoms. 

Alternatively, the produced crystalline Rydberg atom 
configurations can be mapped onto correlated states of 
light and later probed via photon correlation measure- 
ments. Here, a strong 7r-pulse is used to coherently trans- 
fer the Rydberg lattice to a crystal state of atoms in an- 
other hyperfine groundstate \g'). A second laser with 
Rabi frequency Q'(t) drives the transition between \g') 
and a low-lying excited state |e'), which is coupled to \g) 
by a quantum field E(r, t). Analogous to light storage 
schemes in an EIT medium [30] the initially prepared 
quantum state of atomic excitations can be transfered to 
the radiation field E by adiabatically turning on the con- 
trol field Q! . Fig. [5] shows the resulting transfer dynamics 
for the specific example of a three-atom chain, embed- 
ded in a larger condensate [33[. The collective nature 
of the Rydberg excitations not only provides enhanced 
coupling to light but also leads to highly directed photon 
emission along the total wave vector of the applied laser 
fields ktot (fig- 130- The Rydberg chain is mapped onto 
a clear triple-peak structure, with well resolvable delay 
times of ~ lfis. The corresponding second order corre- 
lations, shown in fig. [SJd, clearly reveal the single-photon 
nature of each pulse, as characterized by the absence of 
equal-time photon signals. Since there is a one-to-one 
mapping between the photon density p p h and the initial 
Rydberg state density operator [30[, the depicted time 
correlations directly reflect the strongly correlated nature 
of the prepared many-body state. In this way, correlation 
measurements, provide a direct probe for (im)perfection 
of the underlying Rydberg lattice. 

In summary, we have shown that the excited-state 
number of laser-driven Rydberg gases exhibits a stepwise 
increase as a function of the laser frequency, which bears 
analogies to the Coulomb blockade staircase known from 
nano-scale solid state devices. Our dynamical calcula- 
tions show that frequency chirps allow to adiabatically 
ascend this dipole blockade staircase, thereby forming 
ordered structures of Rydberg atoms out of otherwise 
disordered ensembles. Such crystals naturally realize a 
network of localized collective excitations without the 
need of extended micro- trap arrays. While our simu- 
lations are limited to comparably small excitation num- 
bers, the obtained scaling laws suggest that large crystals 
can be prepared under typical experimental conditions. 
Chirped excitation of lattice-confined atoms to large Ry- 
dberg atom fractions would provide a promising experi- 
mental scenario to study magnetism of long-range inter- 
acting spins and explore their quench dynamics through 
quantum phase transitions. We described how quantum 
optical imaging technique could be used as a real-space 
probe of spatial correlations between the effective spins, 
thereby providing a novel source of non-classical light. 
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